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A new class of braneworld models displaying late-time phantom acceleration without resorting 
to a phantom fluid is presented. In this scenario expansion is fuelled by dark matter together with 
some effective dark energy capable of crossing the phantom divide. Unlike a previous proposal 
of this nature, in these models the effective phantom behaviour remains valid at all redshifts for 
some choices of the free parameters of the models. The construction is based on the generalised 
Chaplygin gas, and the cosmological history interpolates between a standard CDM-like behaviour 
at early times and a de Sitter- like behaviour at late times, so no future singularity is reached. 



: 

Q I. INTRODUCTION 

o : 

Few things in Physics these days can be adjectivated as so intriguing as the reason for cosmological accelerated 
£h ' expansion on large scales. Observations of distant supernovae inform us of a departure from linearity in Hubble's law 
^ for which the most compelling explanation seems to be that the Universe is speeding up. 

The simplest explanation is the presence of a cosmological constant or an exotic fluid making more than two 
thirds of the total matter/energy content of the Universe but there is also the possibility the acceleration is a 
, consequence of infrared modifications of gravity with respect to Einstein's theoretical framework like in the induced 

! , gravity braneworld scenario 0, d, 0, H, @, 0, M, HI, E3, EH, E2, EH EH, EH EB| , which contains as a particular case the 

well-known DGP model 0]. At present we are far from being able to discard either interpretation, so in parallel to 
CIh the observational developments required to answer the cosmic acceleration puzzle, considerable work is being done 
q | along the lines of theory also. 

$_i ■ As the mathematical side of the problem is concerned, finding out about the kinematical/dynamical properties of a 
"qq ' given model for the large scale Universe model starting off from a reconstruction of H(z) is a well posed and clear-cut 
problem. In practice, however, the situation is tougher, as one only has at hand discrete datasets with disturbances 
(for various redshifts) from which one aims at extracting information applicable to all redshift values in a range of 
, interest. Thus, in order to make progress it is vital to smooth out the quantities involved in the description. 
J> ■ Ideally, one would like to do this coarsening using a theoretically grounded description of the background geometry. 
] Given that at present there is not a single model for dark energy which is completely satisfactory, the community 
■ keeps showing interest on new models which firstly do not depart substantially from the characteristics of other dark 
energy models in the market, secondly are likely to give a better fit to the data, and are thirdly are less problematic 
^] . in what theoretical aspects are concerned. 

Observational data do not seem incompatible with dark energy with phantom-like behaviour at z = [l7j . On the 
other hand, this radical behaviour may well be a property acquired by dark energy just recently, so models which do 
the so called phantom divide crossing keep gathering attention (for an overview of the problem see Ref. [l8| ). Quite 
a few realisations of this transition have been discussed in the literature. One can first distinguish general rclativistic 
models, and, among, them there are models with an hybrid of quin tessence and phantom dark energy called quintom 
fl9l . models driven by a fluid called hcsscncc models pll. l22j] which originated after giving quintom models a 
turn of the screw, scalar field models with internal degrees of freedom [23[, k-essence cosmologies (24|, Chaplygin gas 
cosmologies [25| and models sourced by vector- like dark-energ y |26| . Possibilities beyond general relativity are scalar 
tensor theories [27|, f(R) theories [28[, Gauss-Bonnet gravity |29l|. and, finally, explanations which rely somehow on 
the assumption of the existence of extra dimensions (see for example 0, @, 0, H, 0,113, EH, E3 ) ■ 

A set of proposals to achieve that crossing are those made in the context of braneworld models and, in particular, 
those inspired in the LDGP model @,@|- Basically, those authors considered a 4-dimensional brane universe containing 
matter and a cosmological constant on the brane but no phantom fluid [H, Q . The phantom-like behaviour was an 
effect due to the extra 5th dimension. A nice feature of the models is the effective screening of the cosmological 



o 



X 



a Electronic address: 
b Electronic address: 



mariam.bouhmadi@fisica.ist.utl.pt 
ruth.lazkoz@chu.es 



constant, which results in the possibility that its value can be larger that the one observationally allowed in the 
LCDM case. 

The LDGP cosmologies [H, Q where generalised by replacing the cosmological constant with quiessence; i.e. a fluid 
with constant equation of state. These new models are referred to as QDGP cosmologies Q. As a consequence of 
their more dynamical character, the phantom divide crossing is possible, unlike in the LDGP models 1 . This rests 
on the possibility of describing the evolution as driven by an effective dark energy density which together with dust 
is considered to be the fuel of the geometry. Yet, these models can be criticised on the ground that the effective 
phantom picture always breaks down at some point as the effective dark energy density will cease to be positive at a 
sufficiently high redshifts (no matter the values of the free parameters of the model) . 

Here we put forward a possible remedy in the form of an alternative generalisation of the LDGP models, i.e. we 
present a class of models in which the phantom picture stays valid at all redshifts (given some conditions that can be 
met easily). This is achieved by taking the dark energy component on the brane to be a Chaplygin gas (3(j, so we 
can call these cosmologies Chaplygin DGP (hereafter CDGP) models. 

Before proceeding, we recall that the Chaplygin gas by itself can behave as a phantom fluid p3l[ in the standard 
framework of Einstein gravity. However, this phantom fluid could induce a future singularity [32l |. We will show that 
CDGP models also circumvent the future singularity issue that could arise in phantom Chaplygin gas models. Indeed, 
as we will show, in our case the brane is asymptotically dc Sitter. 

In order to show how our models do indeed display the features anticipated in this section, we will provide in 
the following sections analytical details on our construction, and for the sake of completeness, we will also present 
numerical examples to illustrate the effect of changing the values of the free parameters of the models. 



II. THE MODEL 



We consider a DGP_ brancworld model 2 where the dark energy component on the brane is modelled by a generalised 
Chaplygin gas (GCG); i.e. p c h = — A/p" h , and with an extra CDM component. In induced gravity braneworld models, 
which contain as a particular case the DGP scenario, if there is no interaction between the brane and the bulk the 
energy momentum tensor in conserved on the brane, so the usual conservation for the total matter/energy budget 
follows: 

p + 3H(p + p) = 0, (1) 
but restricting it further so that conservation holds for the two components separately we have for the CDM bit 

p m + 3H Pm = 0, (2) 
which integrates to the usual p m = p m o(l + z) 3 , whereas for the dark energy bit we have 

p ch + 3H(p ch + p ch ) = 0, (3) 

which results in 

Pch = Pcw[A s + (1 - A s )(l + z) 3 ^}^. (4) 

The latter is the expression that supports the interest on Chaplygin cosmologies, as it reflects the fact that the energy 
density of such fluids interpolates between dust and a cosmological constant [30l j . As customary the subscript 
corresponds to the current value of a given quantity, and, in addition we define A * = A/(pl^"). 
Putting it all together, the modified Friedmann equation on the brane reads pi, [LH flq 

H/H = vV c + n m (l + z)3 + fU [A s + (1 - A s )(l + 2)3(«+D] ^ - v ^7, (5) 



1 Other generalisations of the LDGP model aiming at describing a mimicry of the phantom divide crossing in induced gravity brancworld 
models can be found in Refs. [ToLUlLll^. 

2 The DGP cosmology has two different brane solutions, corresponding to the two possible ways of embedding the brane in the bulk 

3, 16], the self-accelerating solution (or DGP+ branch) and the DGP- branch. The latter is the one we are considering in the present 
work. 



where 

_ Pm _ gchO o — _ (f\\ 

m - 3H 2 ' ch ~ 3H 2 ' rc ~ Ar 2 H 2 ' W 

are the convenient usual dimcnsionless parameters. 

Throughout the paper we will be assuming < A s < 1 and 1 + a > 0, which imply that the braneworld goes over 
at low rcdshifts to a dc Sitter cosmology (see Eq. ©). 

Similarly to what happened with vacuum energy in the precursor LDGP model 0, we also have a screening 
effect here; the fractional amount of energy in the form of a Chaplygin gas in this CDGP model can be larger than 
in its general relativistic counterpart (which would be obtained on the £l rc — > limit). 

A consequence of Eq. ([5]), which leads to relevant insight on the model, is the constraint 

fim + fich = 1 + Zy/tt^- (7) 

It implies, on the first place that the region £l m + f2 c h < 1 is unphysical (theoretical cuts of this sort are actually quite 
a source of complication when one tries to impose observational constraints [1,0]). Secondly, although the brane is 
spatially flat, the constraint ([7]) can be interpreted in the sense there is a mimicry of closed FRW universes in the 
(f2 m , fl c h) plane. 

Now, a relevant question in this context is whether the brane super-inflates, and in order to answer to it the 
Raychaudhuri equation is of use: 

j=-^^K+ - , (1 ' AsRh(1+ : )3a ^ V (8) 



# 2(H + H ^n^){ [^ + (l-A s )(l + z) 3 ( 1 +«)]■ 

At this point it is in order to stress that the Raychaudhuri equation in braneworld models is modified with respect to 
standard relativistic models, in such a way that the modification of both this equation and the Friedmann equation 
lead to a conservation equation for the matter/energy on the brane which in the absence of interaction between the 
bulk and the brane coincides with the conservation equation one would have in the relativistic situation. 

Recently, an alternative route to the derivation of a cosmological model combining DGP gravity and a Chaplygin 
gas was taken by (33[. In this paper, the DGP-like modification of the Friedmann equation was combined with the 
relativistic form of the Raychaudhuri equation to derive a conservation equation for the possible fluids confined on 
the brane, the equation of state of them being the remaining necessary input. This procedure does not correspond 
to a construction derived from the DGP action (there is though the possibility an action might exist from which this 
derivation could be possible, but as this point has not been addressed in [331 ] . such model can only be regarded as 
phenomcnological at present). 

In these models H is never positive (cf. Eq. ||SJ)), which means the brane never super-accelerates. Note, as well, 
that H vanishes when z — > — 1 while H is positive at this limit, thus reflecting their late-time dc Sitter character. In 
contrast, the large z regime of the models is characterised by 2H = —3H 2 (this can be deduced by taking the z — > oo 
limit in Eqs. © and ©) . As the deceleration factor is well known to be q = — (H/H 2 + 1), and our brane cosmology 
transits between CDM (q = 1/2) and de Sitter (q = —1) limits there is always a certain redshift below which the 
brane accelerates (this redshift is signalled by H = —H 2 ). In Fig. [IJ the dimensionless deceleration parameter q is 
plotted against the redshift z for a fixed set of the parameters Q rc and a. More precisely, we show the evolution of q 
in the left plot of Fig. [T] for a fixed amount of CDM and different values of A s . This plot indicates that the late-time 
acceleration starts sooner the larger is the parameter A s . On the other hand, we show the evolution of q for a fixed 
A s and different amount of CDM in the right plot of Fig. [TJ According to this graph, the larger is f2 m , the later the 
beginning of the cosmic speed up. 

By considering that the Universe is currently accelerating, i.e. H + H 2 > 0, and imposing the constraint given in 
Eq. (JT} it turns out that 

2(n m + n ch ) < i + 3A s n ch . (9) 

This implies the additional constraint 

< A s (10) 



must be satisfied as well. Notice that the last inequality implies that the Chaplygin gas in this model cannot behave 
exactly as dust, i.e. A s ^ 0. 

In the next section we will discuss the main point of the paper, which is the possibility of having not just effective 
phantom behaviour, but also phantom divide crossing. The fact this behaviour is effective, i.e. it is achieved without 
having to resort to scalar fields with a wrong sign in the kinetic energy, allows bypassing the theoretical pathologies 
these models may display at the quantum level. 



Figure 1: Representations of the dimensionless deceleration parameter q against the redshift z. The pair of parameters (n ro ,a) 
has been fixed to (0.01,0.9). In the left plot £l m = 0.2, and the pink solid line, blue dotted line and red dashed-dotted line 
correspond to A s — 0.991, 0.981, 0.971, respectively. In the right plot A s = 0.991, and the pink solid line, blue dotted line and 
red dashed-dotted line correspond to Q m = 0.3,0.25,0.21, respectively. 

III. CROSSING THE PHANTOM DIVIDE 

As said in the introduction, as well as just above, the phantom behaviour in these CDGP cosmologies is an effective 
one, and it rests on the definition of an effective energy density and and effective equation of state parameter w e s, so 
that: 

3H 2 = p m + p eS , (11) 

TJ 

Poff = Pch - 3—, (12) 

l + w eS = --§^. (13) 

It can be shown that the effective energy density does not vanish and is always positive if 

/AO o \ 1+Q 

4A.(1 - A.) > (^F^J • (14) 

Specifically, the latter condition arises from imposing p c g = and the fulfilment of the Fricdmann equation, which 
translates into a relation between p m and p c h, which in turn is equivalent to the vanishing of a quadratic polynomial 
on (1 + z) 3 ( 1+a ' ) . By requiring such equation has no roots, one just obtains expression (|14[) . If this condition is met, 
the sign of 1 + u> c ff is, therefore, completely determined by the sign of p c s or equivalently by 

Pes = Pch - 6^/Ct^H H. (15) 

A straightforwards calculation shows 

9HH*n ch (l-A s )(l + z)^+«) 

Pch = —~ , 5 ■ (10) 

[A s + (1-A s )(l + ^ 1+Q )] 1 + Q 

By combining Eqs. @, ([5]) and (fTB|) . one gets 

PtiS= (H + H 0y /n^) x (17) 

where 

X = il m HoV^-- Hi 1 - A s )n ch (l + zf a [As + (1 - A.)(l + z) 3 ( 1+Q )] . (18) 
At high redshifts (z ;§> 1) the quantity p c g is non-positive as 



X ~ -#(,(! - A s ) — Q ch (l + z)2Jn m + fi ch (l - A s ) — . (19) 
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Figure 2: Contour plots of z c as a function of a and A s for Q ra = 0.04 and Q m = 0.20, 0.25, 0.30, 0.35 (left to right). 



Reversely, when z « — 1 the quantity p ft is non- negative, as one has 

X ~ Q m Ho\/fl r for a > 0. 



(20) 



Thus, by noticing the opposite signs of X at the two rcdshift asymptotic regimes provided fi m ^ 0, it follows that 
this condition is sufficient for p e g to vanish at least once, that is, one is ensuring the realisation of a phantom divide 
crossing. In what follows we will show TO ^ is not only sufficient but also necessary for the passage from the 
conventional to the phantom regime. 

In the absence of CDM (0 TO = 0) the Raychaudhuri equation reads 



2£T = -(1 + w cS )p eS , 



(21) 



and since (in all cases) H < then 1 + w c g > and the phantom divide crossing cannot occur. Just for the record, 
when a Chaplygin gas is advocated to account for both dust and a cosmological constant one speaks about a unified 
dark matter and dark energy scenario, whereas when the Chaplygin gas is considered in the simultaneous presence 
of additional dark matter the Chaplygin gas is regarded as a dark energy component. Our findings demonstrate 
that the Chaplygin gas cannot be the sole component of a DGP braneworld if it is to depict the transition from a 
non-phantom to a phantom-like regime. This seems to indicate that even though the original inspiration for the use 
of the Chaplygin gas was the unification of dark matter and dark energy, applications beyond that original motivation 
are indeed interesting. 

Thus, summarising up to here, we can conclude that the phantom divide crossing is realisable in this scenario in the 
customary effective description, and we want to stress the fact this effective description is well defined at any rcdshift 
if condition (fl"4"|) holds and a is positive. 

It only remains to investigate whether with values of f2 m close to observationally preferred ones, and values of fi rc 
inspired by the closed LCDM "mimicry" (f2 rc <C f2 m ) the crossing will occur at observationally accessible values (say 
z < 1 to be a little strict). This is a tricky task, as one would has on the one hand the difficulty of having to play 
with four parameters (Q m .Q rc , A s , a) and on the other hand the manifest involvement of the pertinent equations. 
Thus, it does not seem possible to give taxative answers in an analytical way. For that reason, we have resorted 
to numerical means and, in particular we have done some contour plots to illustrate the effect of changing those 
parameters. Specifically ,we select constant values of f2 m and Q rc and then draw contours of constant z c (with z c 
denoting the redshift value at which the crossing occurs). For the whole range of parameters considered the conditions 
that the universe is accelerating at present and that p e ff is positive arc satisfied. 

Our choice of range for a is arbitrary and our only reason for the particular values considered is that the examples 
stay close to the original Chaplygin gas (a = 1). In contrast, our choice of range for A s is based on some preliminary 
tests so that our joint choice of values of all four free parameters results in z c values in the close past. Fig. [2] suggests 
a pattern according to which for a fixed value of Q rc typically z c grows with increasing S7 m , whereas, similarly, in the 
view of Fig. [3] one would be inclined to say for a fixed value of f2 TO the typical z c values grows with increasing f2 rc . 




IV. CONCLUSIONS 



Late-time cosmic acceleration is one of the major issues in Cosmology these days and we have addressed it from the 
perspective it is an cxtradimcnsional effect. Wc focus on a particular aspect of it, which has aroused a not disdainable 
activity: the possibility that the equation of state of dark energy (or its effective realisation in our case) has crossed 
the phantom barrier. 

Here, we present an extra dimensions inspired model that is built on the DGP braneworld scenario Q, which 
pioneered the attempts of explaining cosmic acceleration initiated at recent times as an infrared modification of 
gravity 0, 0]. In broad terms the LDGP @, Q and QDGP Q variants of the original DGP_ branch model (see 
footnote [5]) can be viewed effectively as cosmological models filled with CDM and a phantom fluid. The QDGP model 
presents the advantage with respect to the LDGP model that the phantom divide crossing can be realised (so the 
effective fluid can be initially conventional as opposed to its latter phantom- like behaviour). However, both models 
have a drawback: the phantom picture breaks down at sufficiently high redshifts. 

The content of the brane in the LDGP model and the QDGP models is respectively a cosmological constant and 
quiessence, adding in both cases a CDM component to form the whole matter/energy budget. Here we consider a 
model with a Chaplygin gas on top of CDM, and call it CDGP cosmological model. The advantages of this new 
scenario are that (i) the phantom effective picture does not cease to be valid at high redshifts and (ii) a mimicry of 
the phantom divide crossing is also possible, thus saving some of the problems associated to the LDGP and QDGP 
models. We also comment briefly on the consistency problems of another Chaplygin DGP model in the literature [HI]. 

In the preceding sections, we have presented our model starting off from the Friedmann equation and we have 
analysed its main evolutionary features. Specifically, we have investigated the conditions for acceleration to occur, 
and we have concluded that it proceeds provided some restrictions are satisfied by the free parameters of the model. 
We have also addressed the question of super-acceleration, which we show to be unrealisable as in the LDGP and 
QDGP models 3 . Another aspect of the problem to which wc devote considerable attention is the role of CDM, as 
it results crucial for the ability of the model to cross the phantom divide. Finally, we give a numerical account of 
the redshift values at which the crossing occurs, and with the support of some plots, we show there are values of the 
parameters for which the passage from the conventional to the phantom-like behaviour can occur at visible redshifts 
(see Fig. [5] and [3]) while attaining at the same time acceleration at current times (see Fig. [T|). 
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